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Abstract 



We analytically and numerically study the effect of perturbative gluons emission on 
the "Filtering analysis" , which is part of a subjet analysis procedure proposed two years 
ago to possibly identify a low- mass Higgs boson decaying into bb at the LHC. This leads 
us to examine the non-global structure of the resulting perturbative series in the leading 
single-log large- -/Vc approximation, including all-orders numerical results, simple analyt- 
ical approximations to them and comments on the structure of their series expansion. 
We then use these results to semi-analytically optimize the parameters of the Filtering 
analysis so as to suppress as much as possible the effect of underlying event and pile-up 
on the Higgs mass peak reconstruction while keeping the major part of the perturbative 
radiation from the bb dipole. 

1 Introduction 

In recent years, there has been a growing interest in jets studies in order to identify a boosted 
massive particle decaying hadronically, for instance the W boson [H El [3l H], top quarks 
O O [3 [8] , supersymmetric particles [9l [10] and heavy resonances [HI [12l [13] (see also [11] 
for related work on general massive jets). Some of these studies revealed themselves to be 
successful in looking for a boosted light Higgs boson decaying into bb at the LHC [8l[T5 l [T6 | [T7]. 
That of [15[ [16] can be briefly summed up as follows: after having clustered the event with 
a radius R large enough to catch the b and b from the Higgs decay into a sing le jetQ this jet 
can be analysed in 2 steps: 

• A Mass Drop (MD) analysis that allows one to identify the splitting responsible for the 
large jet mass, i.e. separate the b and b and thus measuring the angular distance Rbb 
between them, while suppressing as much QCD background as possible. 

• A Filtering analysis where one reclusters the 2 resulting sub jets with a smaller radius 
and takes the 3 highest-pj subjet^ obtained in order to keep the major part of the 
perturbative radiation while getting rid of as many underlying event (UE) and pile-up 
(PU) particles as possible (used also in [H [TTl [18], and a variant is proposed in [19]). 



^The value chosen was R = 1.2. 

^The value of 3 was found to work well in [15| . 
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Concerning the MD analysis, the only thing we need to know is that we end up with 2 6- 
tagged jets, each with a radius roughly equal to Rbh. Notice that due to angular ordering 
[20l EH [22l [231 ED) these 2 jets should capture the major part of the perturbative radiation 
from the bb dipole. The whole procedure is depicted in figure [1] (taken from |15j). 




Figure 1: The Mass Drop and Filtering analysis in the procedure used to enhance the signal 
from a light Higgs decaying into bb at the LHC. 

In this paper, we are going to focus on the Filtering analysis. One can generalize it with 
respect to its original definition using 2 parameters, Ufut and rjfn^ (as discussed also in [IS]): 
after the MD analysis was carried out, one reclusters the 2 resulting sub jets with a radius 
-Rflit = VfiitRbb and takes the rifut hardest jets obtained. Obviously, the larger the value of 
T/flit the more perturbative radiation we keep, but also the more the UE/PU degrades the 
Higgs peak. The same holds for nf^. So there is a compromise to make between losing too 
much perturbative radiation and being contaminated by soft particles from UE/PU. In |15j . 
the values that were found to give nice results were riflu = 3 and r/AH = min(0.3/-R6;,, 1/2). 
However, these values had been chosen on brief Monte-Carlo event generator study and one 
would like to gain a little more analytical control over them. One question would be for 
instance to understand even approximately how the optimal (nfiit,r/fiit) values change when 
one increases the Higgs ptjj cut, or when the PU becomes more and more important during the 
high luminosity running of the LHC. Though the MD and Filtering analysis were originally 
designed to identify a light Higgs boson, one should be aware that similar calculations may 
apply in other uses of filtering, for instance to study any boosted colorless resonance decaying 
hadronically, including W and Z bosons. 

The article will be devoted in large part to the study of the dependence of the perturbative 
radiation loss with respect to the filtering parameters. As usual in this kind of work, large 
logarithms arise due to soft or collinear gluon emissions, and one is forced to deal with them 
in order to obtain reliable results in the region where the observable is sensitive to this kind of 
emission. We will thus compute analytically the two first orders in the leading soft logarithmic 
(LL) approximation when Uf^i^ = 2 and to all-orders in the large- A'c limilH when n^i^ = 2 or 
3 for small enough values of 7]^^ (section [2]). With these in hand, and using a program 
that allows one to make all-orders leading-log calculations in the large-A'^c approximation, 
we check the analytical results and examine if the small- r/fm limit and/or the truncation of 
the LL expansion can be trusted to estimate the loss of perturbative radiation in practice 
(section [3|) . Finally, in section [H we will analyse the Higgs mass peak width due respectively 
to the loss of perturbative radiation and to the presence of UE/PU, before combining them 
in a simple and approximate but physically reasonable way in order to be able to conclude 
about the optimal parameters choices. 

^Nc — 3 denotes as usual the number of QCD colors. 
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2 Non-Global structure: analytical insights 



2.1 The filtered Higgs mass: a Non-Global observable 

It is now very well known (SgESlETlIMlESlEOlEIlEaESlISlE^ that soft 

or collinear gluons can give rise, in multiscale problems, to the appearance of large logarithms 
in the perturbative expansion of an observable, and more precisely in a region of phase space 
where it is sensitive to the soft or collinear divergences of QCD. In this article, the observable 
considered is AM = Mh — Maitered jet , where Mfiitered jet is the reconstructed Higgs-jet mass 
and Mh is its true mass. AM has the property that it is when no gluon is emitted. We 
are interested in S(AM), the probability for the difference between the reconstructed and 
true Higgs masses to be less than a given AM. In this case, large soft logarithms have to be 
resummed at all-orders to obtain a reliable description of the small AM distribution. 

For this observable, soft gluons emissions lead to powers of In-^^, whereas collinear 
gluons emissions leads to powers of In-^^. In this study, gluons are strongly ordered in 
energy (the first emitted gluon being the most energetic one, and so on), and we aim to 

control the ( Ug In ) series, in a region where 



lni^»ln:^. (1) 

AM iifl,, ^ ^ 

Therefore, at leading-log accuracy, one has to resum terms like 

where all the fk are functions to be computed. We thus disregard all the subleading terms, i.e. 
those suppressed by at least one power of In Unfortunately, such a calculation is highly 
non-trivial due to the fact that the observable is non-global. This property, first studied in 
|35j . means that it is sensitive to radiation in only a part of the phase space. In the case of 
AM, only emissions of gluons outside the filtered jets region contribute to the observable (cf 
figure [2]). As a consequence of this property, one must consider soft gluons emissions not just 
from the bb dipole (usually called primary emissions, the only ones that would be present in 
QED) but also from the whole ensemble of already emitted gluons [351141]. As the number of 
gluons is increased, the geometry and the color structure of all these gluons become rapidly 
too complex to perform an analytical calculation. Therefore, to deal with this, one is forced 
to apply numerical Monte-Carlo calculations that can easily take care of the geometry. But 
the colour structure remains prohibitive, and one must usually also resort to the large-A^c 
approximation in order to go beyond the 2 first orders in perturbation theory \35\ [411 [42l [43] 
(though some authors have derived some analytical results in special cases [44 [ US] and others 
have examined contributions beyond the leading large- A'^c approximation |46[ [47] ) . 

However, before considering numerical calculations, some results can be derived analyti- 
cally at 2'^'^ order for = 2 (where /i and /2 are computed exactly) and nan = 3 (where 
only the leading behaviour of the fk in In and Nc is looked for). 
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(a) 




(b) 



Figure 2: Configurations leading to non-global logarithms when (a) Ufut = 2 and |(b)| nfiit = 3. 
In each case, the hardest gluon 1, which is inside the filtered jet region, emits a softer gluon 
2 outside the filtered jet region. 

2.2 Some results for n^^ = 2 

Perturb at ively, one can write I](AM) as 



S(AM) = 1 + ^4(AM) 



(3) 



k=l 



where /^(AM) is the O (a^) contribution to the observable. To simplify the calculation, 
S(AM) will be computed using the anti-Zcj algorithm [l8], even if the numerical study will be 
done using the C/A algorithm [491 150| to be in accordance with the choice in |15j. However, 
the anti-kt algorithm is enough to catch the dominant behaviour of the leading-log series, in 
the sense that it does not affect the leading large collinear logarithm in the function at 
small Rfiit- 



fk 



Rbb 



ak In'^' 



Rbb 



+ O i In'^-i 



Rbb 



(4) 



, Rmt ) \ Rftit / V \ ^flit , 

i.e. Ofc is unchanged when moving from C j A to anti-A:j0 This jet algorithm gives simpler 
results because the gluons outside the filtered jet region tend not to cluster with the ones 
inside. It is this property which ensures that the hardest jets in an event are generally perfect 
cones, as particles usually cluster with the hardest ones in their neighbourhood first |48j . 

As a first step, primary emissions are considered, defined to be those one would obtain if 
gluons were only emitted from the hb dipole (as for photons in QED). 



2.2.1 Primary emissions 

Due to the use of the anti-fcj algorithm, the result of the primary emissions can be easily 
shown to exponentiate, as will be roughly seen in the next section with the O (a^) analysis. 

''When i?fiit ~ the discarding of the O ^In*""^ t'^™^ '^"^ "^^^ ^ priori justified, but fig. O 

which compares numerical results obtained using C jA with analytical estimates using anti-fct, supports the 
dominance of the leading collinear logarithms. 
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Here, we just review the very well known result that the contribution to S(AM) from primary 
emissions, denoted S(^)(AM), can be written as| 



S(^)(AM) = e^i(^*^), (5) 



with: 



ii(AM) = J -J-h_Mih) (e (k, G Jtb) + e (fci ^ e (am - AM{h)) - 1 

(6) 

M{ki) is the matrix element squared for emitting one soft gluon from the bb dipole (the b 
quark is taken to be massless): 



Concerning the notations, ^A;i G J^^ j equals 1 when gluon 1 is emitted inside the jet regions 
around b and b, denoted by J^g (and is otherwise), which, for Rait < Rbb, is just 2 cones of 



radius Rfat centered on b and b (figure 2(a) ). Then, concerning the expression in brackets in 
eq. ([6]), we separate the 2 different regions where the gluon can be: either inside or outside 
the filtered Higgs jet. The first term Q{ki £ Jj^i) means that the gluon does not contribute to 
the observable (as it is kept in the Higgs jet, the reconstructed Higgs mass is the true Higgs 
mass: AM{ki) = 0). If the gluon is outside the filtered jet region (second term), then it does 
contribute to the observable: 

AM(fc) ~ fc^— ^ , (8) 

PtH 

up to prefactors that can be neglected in the leading-log approximation, see appendix lAl 
Finally, the —1 stands for the virtual corrections, for which there's obviously no loss of mass 
for the Higgs, and whose matrix element is just the opposite of the soft real onel§ One thus 
obtains: 

Ii{AM) = - ^M(fci)e {AM{ki) - AM) . (9) 

Jki(J,i (27r)32|A;i| 

The computation of this integral in the boosted regime, where ptn ^ Mh, or equivalently 
Rbb <C 1, is done in appendix[Al From now on, we will essentially use t/hh = Rfn^/Rbb instead 
of i?fiit and we define n = n^n and r] = for more clarity in mathematical formulae. In order 
to keep in mind that it depends on the 2 parameters of the Filtering analysis, the distribution 
S(AM) is renamed T^'^^^i], AM). What we obtain at fixed coupling is the followinglll 

AM) = e-^^^")!-^ , (10) 



®The superfix (P) serves as a reminder that only primary emissions are being accounted for. 

^Even if the result seems obvious here, this way of doing the calculation can be easily generalised to higher 
orders and other kinds of jet algorithms. 

^To obtain the result at running coupling, one simply makes the replacement (see eqs. (|29I30[) later in the 
article) : 

, Mh I , f 1 
a. In . . , — >■ — — m 



AM 2/3o \l^2l3oasiMH)ln§§ 



5 



with 



J(ry) =2 In (i-f 



TT 



— 1) 



/li-l /2u+l\ 1 
arctan I t\I r I ii - < r? < 1 . 



u + l\ 2u-l 



fill 



(12) 



We give the value of J(l), a quantity that is important to discuss some aspects of the results 
obtained in the following sections: 

J(l) - 0.646. (13) 

Notice that the case r/ > 1 will not be used in this study, but is mentioned in appendix lAl 
The function J{r]) is plotted in figure [3l 




0.4 0.6 0.8 1 

Figure 3: The coefficient J{r]) in front of the primary soft logarithm In 
Two remarks can be made: 

1. The result does not depend on the energy fraction z of the Higgs splitting into bb. 

2. When r/ <C 1, another large logarithm In i appears due to collinear enhancement. 

2.2.2 Non-global contributions 

Now, we turn to the O (a^) term, and more precisely to the contribution of the non-global 
terms that have to be added to the primary logarithms computed in the previous section. 
That corresponds to the analysis of /2(AM) in the perturbative expansion of E(AM) from 
eq. dSD . The matrix element squared for 2 real gluons emission is expressed as [Ml ESI ED E2] : 



M{ki real, k2 real) = (AnasfiWi + W2) , 



with Wi = ACl 



(Pb-Pb) 



(Pb-Pb) 



W2 = 2CfCa 



{pb-ki)iki-Pb) {pb-k2){k2-Pb) ' 
iPb-Pb) ( {Pb-h] 



(14) 
(15) 



+ 



(Pb-ki) 



{Pb-Pb) 



{pb.ki){h-Pb) \{Pb-k2){k2-ki) {Pl-k2){k2.ki) {pb-k2){k2-Pb) 

(16) 
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This expression is valid when there is a strong energy ordering between the two real gluons 
1 and 2, either Ei ^ E2 or E2 ^ Ei (the formula is completely symmetric under the 
interchange ki -f-T- k2)- For the cases with one or both gluons being virtual, the following 
matrix elements are obtained, valid only when Ei S> E2 [52] : 

M{ki real, /c2 virt) = -{AnasfiWi + W2) , 
M{ki virt,/c2 real) = -{4Tras)'^Wi , 

M{ki virt, /C2 virt) = {inasfWi . (17) 

Using these properties, separating the 4 phase space regions depending on whether the gluons 
are inside or outside the filtered jet region in the same way as was done for Ii, and defining 
dk as: 

dk = (18) 
(27r)32|A;| 

we can then write I2 in the following form: 

dkidk2{4^asfQ{Ei-E2){ 

@{k^ G .hl)@{k2 e %)((TVi + W2) - {Wi + W2) -Wi + Wi) 
+ G(A;i G Jfefe)e(A:2 i Ji,i){{Wi + W^2)e(AM - AM(A:2)) - {Wi + W2) - VFiG(AM - AM(/c2)) + Wi) 
+ G(A;i i Jbb)eik2 G ^fe5)((H^i + W2)e(AAf - AM(A:i)) - (H^i + W2)e(AM - AM{ki)) - Wi + Wi) 
+ G(A;i ^ Jfe5)e(A;2 ^ Ji,5)((VKi + i^2)e(AM - AM(A:i, A;2)) - {Wi + T^2)e(AM - AM{ki)) 

- Wie{AM - AM{k2)) + Wi)] . (19) 

For each phase space region, the 4 terms (/ci,/c2) = (real, real) — (real, virt) — (virt, real) + 
(virt, virt) are considered. The strong energy ordering Ei ^ E2 implies that AM(/ci,/c2) = 
AM(/ci), and one immediately gets: 

h= j dkidk2{A-nasfQ{Ei - ^2)6(^1 ^ Jbi)Q{k2 i .hi)Wi@{AM{k2) - AM) 

dkidk2{4T^asfQ{Ei - E2)Q{ki e Jki)Q{k2 i %)W2G(AM(A:2) - AM) , 
= /^^)(AM) + 4^^^(AM), (20) 

where l'^^ (AM) corresponds to the first integral containing the function Wi whereas l'^^'^ (AM) 
corresponds to the second integral with the function W2. I2 is just the second order contri- 
bution to the primary emissions, already computed above. To be convinced, one can notice 
that {A-KUsfWi can be expressed as the product of 2 one-gluon matrix elements M{ki)M{k2) 
and, when Ei ':$> E2, 

e (AM(A:2) - AM) = 6 {AM{k2) - AM) 6 (AM(A:i) - AM) , (21) 
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(P) 

if ki and k2 belong to the same phase space region. Therefore I2 can be written in a more 
symmetric way: 

iPiAM) = ]^ j dkidk2{AT^asf@{kr i JbbMh ^ Jbb)WieiAM{ki) - AM)G(AM(A:2) - AM) , 
= \(^j dke{k i J^i)M{k)e{AM{k) - AM)^ , 

= i(Il(AM))^ (22) 

so that it corresponds to the second order perturbative expansion of the result eq. (0), ob- 
tained with primary emissions only. 

The important term for this section is the one containing W2, denoted by As 
mentioned in section 12.11 it receives a non-zero contribution when the hardest gluon 1 is 
emitted inside the filtered jet region whereas the softest gluon 2 is emitted outside. For the 
opposite configuration, there is an exact cancellation between gluon 2 being real and virtual. 
Here again the computation of ^2^'^^ is postponed to appendix [Al giving directly what will 
help to interpret some results later. S2 is defined such that 

ir\ri,AM) = \C,C, (^In [^^)^ S2{,) , (23) 



where we explicitly introduce the dependence on 77 and we factorize out the soft divergence, 

M„ 
AM- 



still revealed in the large logarithm ln^#. When rj < 1/2, the result for ^2 can be written 



as: 

3 Jo ^^1 Jo 



ui Jo U2 \^{l-r]^{ui + ui)Y -Ar^^uiui l-^^K + ^l)/ ' 
+ 4//^ + Urf + O [rf ] . (24) 



The important point to notice in this result is the absence of collinear logarithms, which 
would appear as In^, contrary to the primary emission case (eq. (jlip ). So that the primary 
emissions dominate for this observable, at least for 77 sufficiently small. 

As mentioned in previous studies |35l I53j . one notices the presence of "vr^ terms" in 
non-global results at second order. 



2.3 Some results for riflit = 3 

The goal in this section is to have an estimate of the analytical behavior in the large Nc limit 
of S(")(r7, AM) for n = 3, which is the probability of having no second gluon emission leading 
to a AM' greater than AM. Notice that, contrary to the previous part where we obtained 
the function S^^^, only the leading behavior in L = In ^ and N^. will be derived, so that in 

this context S*^^-* can be simply writtenH 

s(2)(L,f) = e-^^=^* (25) 

*This results simply from the combination of equations pop and ([TT} with 77 ^ 1, and Cf = in the 
large limit. 



8 



where for further convenience we introduce the parameter t = ^ In and we change the 
arguments of S which becomes now a function of L and t. In this formula, 2L = 2 In 

^ _ ' -Kflit 

can be interpreted as the "logarithmic size" of the bh dipole, i.e. the allowed phase space in 
rapidity for an emission from this dipole (in its center of mass frame) outside the jet region. 
The parameter t means that this emission cannot occur with a t' between and t. 

Now we turn to Yi^^\L^t). To have no second gluon emission in [0,t], either there is no 
first gluon emission in [0,t] outside the jet region (which corresponds to S('^)(L,t)), or there 
is such an emission but the new dipole configuration is prohibited from emitting a second 
gluon in [0, t] outside the jet region. This is depicted in figure HI As the calculation is done in 



s(3)(L,t) ^ . ^ ^ r 



dldt' 




Figure 4: How to compute the leading behavior of S('^)(L, t) from S(^)(L, t) when L S> 1 and 
Nc'^ 1. In the second term, t' is the gluon's emission scale. 

the large- A'^c limit, after the emission of a first gluon, the second one cannot be emitted from 
the bh dipole, but only from the bg and bg ones. Fig. [5] can be translated mathematically as: 



s(3)(L,t)~s(2)(L,t)+ f dt'ANj:^'^\L,t') dlT,^'^\L,t-t')^ 

Jo Jo 



(2) 



il,t-t'). (26) 



Notice that Lbg, the logarithmic size of the bg dipole in figure [H does not depend on / in the 
leading collinear log approximationlfl In this expression, 4A^c^S(2)(L,f)dt' is the probability 
not to emit the first gluon in [0, t'] and to emit it only at t G [t' , t' + dt']. The remaining part 
Jq dl S(^)(L, t — t')E(^)(/, t — t') is the probability to emit no second gluon from the bg and 
bg dipoles in [t',t]. Using eq. (i25]l for S^^^ T.^^^ is then given by: 



/ MNcLt 1 -t' \ 

s(3)(L,t) ^ e-^^^^* f 1 + / dt'^-^ . (27) 







One can easily show the following relation: 



Lbg^2L + 0{ e'-"- 



If we introduce the neglected component of Lbg in the S''^' calculation eq. (|26p. then we have to compute an 
integral of the form 

/■I 

dl- 



l + 0{e''-^) 

Expanding the exponential and keeping the term of order k gives 



The leading O {{Ltf) term is already taken into account in eq. (|27p . Therefore, including the I dependent 
component of L^g gives rise to terms of the form N^L^~'^t^ at order k, suppressed by 2 powers of L with 
respect to the leading one. 
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Two limits can be considered: 

r I - liiNcLtf + O {{mcLtf + Net) [fmcLt<.l, 

(28) 

[ e-"^^-^^ {ln{4NcLt) + 0{l)) if 4iVcLi > 1 . 

The limit ANcLt ^ 1 reveals two important aspects: 

1. One can notice the absence of the 0{Lt) term, which is indeed the goal of the filtering 
analysis as it was presented in its original version |15| : it is intended to catch the major 
part of the 0{as) perturbative radiation. It cannot catch all the 0{as) contribution 
because a hard gluon emitted at an angle 9 > Rhi, from the b and b escapes the filtering 
process as it is rejected by the Mass Drop analysis. Therefore, when ANcLt <C 1, the 
expansion eq. (j28p misses a term 0{Nct), but this is legitimate in a leading coUinear 
log estimate. Notice that the missing term is simply —J{l)Nct where J(??) was given 
in eq. (112]) . 

2. it shows that the purely non-global result for n = 3 contains large collinear logarithms 
L, contrary to the case n = 2 (eg. (1241) ). Indeed, the primary result for n = 3 at second 
order can be proven to behave aq^j —32Cp{Lt)'^ at order a^, so that the 5*2 term for 
n = 3 should be equivalent to — 8Ci?C^(Lt)^ at large L. 

Having understood some analytical features of the Filtering analysis, we now examine 
what can be learnt from a numerical calculation of the reconstructed Higgs mass observable. 



3 Non- Global structure: numerical results 

In all that follows t is defined so as to gather all the information about the soft logarithms 
in a running coupling framework: 

2vr Jo h V PtH J 

1 fP'H dkt / Mh 
a., ( kt 



47r^o^'' (l-2/3oa3(MH)ln^ 



I . 1 ' (29) 



where the last equality holds at the one- loop level and /3o = ~~j2^^~^- The argument of Ug 
was taken as the gluon's transverse momentum with respect to the Higgs boson direction, of 
order kt^^, kt being its transverse momentum with respect to the beam. In the case of a 
fixed coupling constant Og, the definition for t here coincides with that of section [2^31 

^ = -1-^- (30) 
27r AM ^ ^ 



^"in fact, one can show the following general estimate for the primary emissions in the leading soft and 
collinear approximations: 

(L, t) = e-*''^" J2 ^^^^Ir^ ■ 

k = 
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But from now on, and unless stated otherwise, t is given in the running coupling framework, 
eq. ([29|) . and the function T,{r],AM) is rewritten as 

To get an idea of the range of values covered by t, table [D presents a few t values corre- 
sponding to a given AM for a Higgs mass of 115 GeV {as{M}{) = 0.114). It reveals that the 
physical values for t are below 0.15. 



AM (GeV) 


1 


2 


5 


10 


20 


50 


115 


t 


0.141 


0.108 


0.075 


0.054 


0.036 


0.016 






Table 1: Correspondence between AM and t for some particular values. 



To numerically investigate non-global observables, two approaches can be followed: 

• an all-orders approach where one resums the leading-logs at all-orders in the large- 
Nc limit, the output being the function S(t), i.e. the probability that the loss of 
perturbative emission results in a Higgs mass in the range [Mh — AM{t), Mh], with 

AM{t) = MHe 2/3o"»^' (31) 
simply obtained by inverting the relation eq. ()29p . 

• a fixed-order approach where the first few coefficients from the expansion of S(t) are 
computed in the large- A'^c limit. More precisely, if S(t) = X] H i^ct) , then the pro- 

A:=0 

gram returns the first few coefficients Cfc. 

Prom a numerical point of view, the way to write an all-orders program was explained in 
[35) . On the other side, a result at fixed-order may be obtained by developping a systematic 
approach like the one presented at second order in eq. ()19p . For the filtered Higgs jet mass 
observable, we used the Fastjet package [54] to perform the clustering (and mass-drop -|- 
filtering) with the C/A algorithm, consistently with the choice made in |15j . 

As the all-orders program gives immediately what we are looking for, which is S(t), we will 
use it (section |4.1|) to compute the perturbative Higgs width. But in order to check it and be 
confident with the results obtained, we compare them with the previous analytical estimates 
and see how well the perturbative leading log series fits them. This leads us to study the 
behaviour of the higher order terms and to gain a better understanding of the convergence 
and structure of the non-global series. Though treated in more details in appendix [UJ the 
main points are mentioned in this section. 

3.1 Comparison with analytics 

Using the all-orders Monte-Carlo program, a comparison between the all-orders numerical 
curves obtained using the C/A algorithm and their corresponding analytical estimates ob- 
tained previously with anti-fc^ in eqs. (|25|27p can be done. The results are presented in figure[5] 
and show good agreement, at least in the region of physical t values. 

Notice that the slight discrepancy between analytical estimations and numerics starts to 
occur at t > 0.1, which is at the edge of the physical region (cf tabled]), beyond which AM 
would be below the perturbative scale of around 1 GeV. This agreement manifests that: 
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Figure 5: Comparison between numerical all-orders result (obtained using C/A algorithm) 
and leading collinear logarithm estimate of derived with anti-Zct for (a) n 
n = 3 for 2 values of r]: 0.1 and 0.5. 



2 and (b) 



• In the physical region, the leading terms in [us Ltf, with L = Ini seem to completely 
dominate and we do not need to compute the subleading corrections. 

• One can use these analytical expressions to get an accurate estimate of the reconstructed 
Higgs peak width. 



3.2 Comparison with fixed-order results 

The structure of the non-global series at fixed-order is now examined so as to independently 
cross-check the all-orders program and to understand if the perturbative leading-log series 
can be usefully truncated. 

As an example, figure [6] compares the all-orders result to the fixed-order ones up to 
for n = 2 and two different values of rj (only the coefficients with an uncertainty of at most 
a few percent are plotte<£l]). The curves are represented up to t = 0.3, which is far beyond 
the physical region but is instructive to study the convergence of the series. 

The left plot for = 0.3 shows a nice convergence of the perturbative series eq. ([3]), as 
the t range for which the all-orders and fixed-order curves coincide grows with k. However, 
the second plot for r] = 0.9 gives an unexpected result: the fourth order diverges with respect 
to the third one, in the sense that the point of disagreement is shifted to smaller t. The 
question arises whether this divergence will remain at higher orders. To answer it, one needs 
to go further in perturbation theory. In appendix [Cl a parallel is made between the filtered 
Higgs jet observable and the slice observable, studied for instance in 01], for which, due to 
computationnal speed, it is possible to obtain reliable coefficients up to order 6. The same 
effect is observed and is even enhanced at orders 5 and 6. Therefore, it seems that the fixed- 
order information cannot be safely used in general: one has to be aware that the leading-log 
large- A'^E non-global series may be divergent for any value of t. 

^^This uncertainty obviously increases with the perturbative order, but also with 77 because at small r), the 
coefficients are sensitive to the large logarithm In - , which is easy to compute. 
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Figure 6: Comparison between fixed-order (FO) expansion and all-orders result when n = 2 



for (a) r] = 0.3 and |(b)| r] = 0.9. Both fixed-order and all-orders results were obtained using 



the C/A algorithm. 



4 Choice of the filtering parameters 

In the previous sections we examined the structure and convergence of the perturbative 
leading-log series, analytically and numerically. We could then cross-check the analytical 
expressions and the fixed-order approach with the all-orders program, which we are going to 
use throughout this part. 

We would like to decide how one should choose the filtering parameters {n,rj) depending 
on the level of UE and PU as well as the pt of the Higgs boson. Here, we do not claim to make 
an exact and complete analysis, but we want to obtain some estimates. First, we consider the 
width of the Higgs mass distribution separately in presence of perturbative radiation (using 
the all-orders results) and UE/PU (using a simple model for it). Then, we try to minimize 
the Higgs width in presence of both of these effects. Finally, we will estimate hadronisation 
corrections. 

In all this part, we set the Higgs mass Mh at 115 GeV, as in |15] . 



4.1 Study of the Higgs perturbative width 

As we could see in the previous sections, even without considering additional particles from 
UE/PU, AM = Mh — Mflit^cdjot 7^ because of the loss of perturbative radiation. The 
Higgs boson thus acquires a perturbative width, denoted 5MpT- At first sight, knowing the 
distributions i;(")(r/, AM), one might simply define it as: 

SMpT = 2v^(AM2) - (AM)2 , (32) 

as we do for gaussian distributions for instance. Unfortunately, if we simply take n = 2 as an 
example and if we consider the primary emission result eq. (jlOp , we can deduce the following 
distribution for AM: 

dS(2)(r?,AM) ^ asC{r,) 1 

dAM j^^sCiv) ^M^-»sC{rj) ' ^"^"^^ 
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with 



Civ) 



CpJiv) 



vr 



Computing (AM) and (AM^) implies dealing with integrals of the form 

A/h W a l\,f f-Mff ^ ^ 

, ^^^i!L_AM = / dAM AM°=^('') , 
AAfi-^-c-Cr,) 



dAM 



AAf2 



Mh 



(34) 

(35) 
(36) 



Such integrals give a large importance to the AM ~ Mh/2 region, where there should be very 
few events, and do not describe what happens in the neighbourhood of the peak near AM = 0. 
Therefore, the definition eq. (|32p does not seem adequate for the perturbative width. That's 
why we shall adopt another definition, adapted from [55]. The Higgs perturbative width is 
defined as the size SMpx for which a given fraction / of events satisfy < AM < 5MpT- 
Using the all-orders function previously computed, this is equivalent to solving the equation 
S(")(r/, AM) = /. This leads to the width function 6MpT{n, rj, /). 
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Figure 7: Perturbative width of the Higgs boson (in GeV) as a function of r] for several values 
of n when (a) / = 0.68 and|(b)| / 



0.95. 



Fig.[7]shows 6MpT as a function of forn = 2 . . . 6. When AM ~ 50 GeV (i.e. ~ Mh/2), 
one should be aware that soft approximation loses sense and results on these plots should 
no longer be taken seriously. We chose the values / = 0.68 and / = 0.95, corresponding 
respectively to 2a and 4cr for gaussian distributions, to show that the Higgs mass perturbative 
distribution is not gaussian (otherwise, going from 2a to 4o" would have multiplied the width 
by a factor of 2, see also eq. (j33p ). One important thing to notice is a kind of "saturation" 
effect that one observes for rj close enough to 1 for every fraction /. It manifests itself as a 
flat curve at a value 6MpT = 5Msat{f), independent of n. For instance, 5Msat{f = 0.68) ~ 1 
GeV and SMgatif = 0.95) ~ 33 GeV. This can be understood simply by considering that 
when the radius of the filtering is large enough, say r] > r]sat{n)-, it captures (almost) all 
the particles resulting from the Mass Drop analysis, i.e. all those that are within angular 
distance i^^ft from h or 6, but it still fails to capture particles outside the Mass Drop region 



^The probability to emit a gluon outside the MD region in [0, t] is roughly given by 1 — e 



-J{l)Nat 
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Of course, the larger n, the smaher rjsatin) as we keep more jets. This saturation property is 
equivalent to saying that all the functions S(")(r/, AM) become independent of n and t] when 
r] > r]sat{n). 

For the rest of this analysis we keep the value / = 0.68, even if it is not clear which 
value should be chosen, and more generally what should be the relevant definition of the 
Higgs perturbative width. However, we will mention in section what happens if we vary / 
between / = 0.5 and / = 0.8, so as to obtain a measure of the uncertainty of the calculations. 

The curves in figure [7] only give us an overview of the scales involved in the Higgs boson 
width. But one can go a little further. At small r], we should get a large collinear enhancement 
revealing itself as a large logarithm L = In ^ multiplying t. The perturbative expansion is 

thus a series in {NcLtf. As a direct consequence, at small r/, the all-orders function S'^"'-*(r/, t) 
can be written as a function of a single variable Ti'^'^\NcLt). Solving the "width equation" 

EW(iV,Lt) = /, (37) 

gives 

tpT = ^^^, (38) 
where tpT is simply related to 6MpT by 



1 



tpT = — - In j-f— , (39) 

4vr/3o \l-2p,as{MH)\n^) 

and where CpT{n,f) is a function, independent of rj, which increases with n and decreases 



when / increases. This is confirmed by figure 8(a) which shows that tppL is indeed indepen- 
dent of 77 as long as rj and n are not too large. 

As an example, for n = 2, let us take the simple result T,^'^\L,t) = ^-^^^^^ from eq. ()25p 
in the small r/ limit. It was shown in section [3T] that this result is very close to the all-orders 
one in the physical t region. Solving Ti^'^\L,t) = f immediately implies 

C^pr(2,/) = ^, (40) 



which, for / = 0.68, gives Cpt — 0.032 in accordance with figure 8(i 

One observes that tpxL is not strictly speaking a constant for higher n values. This may 
be due to the saturation effects discussed above. Indeed, even at large L, the perturbative 
expansion is not only a function of Lt but also of t for the lowest orders, as mentioned at the 
end of section 12.31 

n—2 +QO 

sW(L,t) = l + ^afct'=+ {ak{Ltf + o{L^-h'^)), (41) 

k=l k=n—l 

If we only had QED like emissions, i.e. primary ones, with the use of the anti-fc^ jet algorithm, 
we would obtain at = ^ '^^fci^''^ for k < n — 2, where J(r/) was derived in section 12.21 As 
n increases, the term ait becomes more and more important with respect to a„_i(Lt)"~^, 
leading to larger and larger deviations from the simple law tp^L = constant. However, until 
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Figure 8: (a) tpxL as a function of t] for / = 0.68 and different values of n. The saturation 
curve simply comes from the fact that all widths saturate to the same constant SMgat for 77 
large enough, and its equation is therefore given by tsatL (with tgatif = 0.68) ~ 0.136). |(b)| 
6MpT as a function of rj for / = 0.68 and different values of n (curves with points). For each 
n is also represented the corresponding approximate width (lines) given by eq. (j42p . 



n 


2 


3 


4 


5 


CpT 


0.032 


0.078 


0.117 


0.149 


Vsat 


0.79 


0.56 


0.42 


0.34 



Table 2: Cpt and ijsat as a function of n when / = 0.68. 



n = 5, assuming IptL is a constant at small r] seems a good approximation. Therefore, using 
eq. (|39|) . one can model the Higgs perturbative width in the following form: 

i_ 



6MpT{n,L,f) 



l-e 



) z 



lirj < risat{n,f) , 
if 7? > rjsatinj) . 



(42) 



Vsatin, f) is given by the intersection between the curve tpT = Cpt/L and tpx = tgat- 
Therefore: 

Vsat{n,f) = e . (43) 

Table [2] shows Cpt and rjsat for / = 0.68 and different n values. Figure |8(b)| shows the 
curves corresponding to the parametrisation eq. (j42p . We can see that it works rather well 
for all values of n except n = 2 in the region rj ~ 0.4 — 0.6. This can be improved using the 
relation J{r])tpT = constant, which works better for n = 2 because it is exact for primary 
emissions with anti-kt- But implementing it would not change the main conclusions presented 
in sections 14. 3114. 51 Therefore, for the sake of simplicity, we will not use it here: we keep 
eq. (^2!) as the expression for SMpt for the rest of this study. 

Of course, were it only for the perturbative radiation, it would be nicer to choose > rjsat 
in order to catch as many gluons as possible, leading to 6MpT — )• SMgat- But we also 
have to take into account Initial State Radiation (ISR) from the incoming qq pair and non- 
perturbative effects like PU and UE that can spoil our Higgs neighbourhood, thus increasing 
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the jet mass. 

For the purpose of this article we wih only add UE and PU to the Final-State Radiation 
(FSR) effect studied above. We will thus ignore ISR, this partially for a question of simplicity 
of the analysis, but also because the results of work such as [181 ES] suggest that for LHC 
processes whose hard scales are few hundred GeV, the crucial interplay is that between FSR 
and UE/PU. This is evident in the preference for small R values in dijet mass reconstructions 
in those references, where ISR is not playing a major role. Similarly we believe that the 
optimal values of rj that we will determine here will have limited impact from ISR, though 
we shall not check this explicitly. 



4.2 Study of the Higgs width due to underlying event and pile-up 

For this simple analysis, which does not aim to give precise numbers but only an estimate of 
the influence of the UE/PU on the mass of the Higgs jet when we vary for instance r], n, or 
when the Higgs boson becomes more and more boosted, we model the UE and PU as soft 
particles uniformly distributed in the {y,(p) plane \57\ I58j . and with transverse momentum 
per unit area denoted by p. In order to get this estimate, we consider the simple case of 
a symmetric (z = 1/2) Higgs decay along the x axis. In the limit Mh <^ Ptu^ the Higgs 
momentum pH is given by: 

PH={ptu + ^.PtH^^^^ ■ (44) 

The UE/PU momentum, denoted pt/sEl is simply the sum of all the UE/PU particles g 
belonging to the filtered jet J. Still in the limit Mh <^ pt^, we recall the following formula 

mi: 

U^^^^ME. ,45) 

Throughout this section, we will apply it with z = 1/2. We can now write AM = Maitered jet — 
Mh as: 

AM=^{(pH+PUEf-M],) , 

In the last line we used the approximation: 

OgH -0,H = ^, (47) 

which comes from the fact that the UE and PU particles tend to cluster around the pertur- 
bative radiation, which is usually close to the b and b because of the collinear logarithmic 



^■^For brevity, we define puE to be the sum of the UE and/or PU particles' momentum but without refer- 
encing the PU dependence, which will always be implicit. 



17 



divergence of QCD. As all the filtered UE/PU particles flow approximately in the same direc- 
tion, the remaining sum is just the total transverse momentum of the UE which, by definition 
of /3, is equal to pA, A being the total area of the filtered jets0 We thus obtain 

AM , (48) 

PtH 

with 

{A) UTTT^^Rl , (49) 

for the C/A jet algorithm, taking into account the anomalous dimension that comes from the 
fact that there should be some perturbative radiation in the jets (cf figure 14 in [57]). Notice 
that eq. (j49p is only true if all the jets do not overlap, so usually when r] is small enough. 
But this is sufficient for the purpose of our study, and we shall use this formula in all the 
following calculations. The correction eq. ()48p for AM only induces a shift towards higher 
masses of the Higgs mass peak. However, there are 3 sources of fiuctuations that give a width 
to this Higgs peak: 

1. p is not strictly uniform in the {y,4>) plane in a given event. 

2. p is not the same from one event to the next. 

3. The jets' area fluctuates. 

Following , we can write the total UE /PU transverse momentum contributing to the Higgs 
Pt as 

pt^^ =pA± (yja + ASp + pE^ , (50) 

where 

^ — \l (P^) ~ (p)^ with (...) a spatial average in a given event , (51) 



bp = \l (/j2) — (p)2 with (...) an average over events , (52) 
S = V(^^> - k^Y with {A) the average over events of the filtered jets' area . (53) 

For pure UE events, i.e. without PU, these terms can be estimated [571 [58]: 

^r/£; ~ 2 - 3 GeV/area , (54) 

a\jE - 0.6puE , (55) 

SpuE =i O.SpuE , (56) 

S ~ 0.26^yn^Tr]'^Rl^ . (57) 

Though puE seems to be around 2 GeV/area, the tuning used in [15] was closer to 3 GeV/area, 
the value that we choose here. In presence of PU, i.e. when there is more than 1 pp collision 
per bunch crossing at the LHC (thus leading to the emission of other soft particles), p, a 
and 5p have to be modified. We define Npu to be the number of pp collisions in a bunch 



*in the active sense, see [57] . 
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crossing except the one at the origin of the hard interaction. We use a simple model to write 
the parameters of the UE/PU as: 



1 + ^) PUE , (58) 
^ - \/ 1 + -^'^UE , (59) 



Sp^^l + ^SpuE. (60) 

Some comments are needed: since p measures the level of noise, it should grow like Npjj. 
In the expression 1 + Npu/A, the 1 corresponds to the pp collision that leads to the UE 
and to the hard interaction, whereas the Npu/A term simply corresponds to the other pp 
interactions and could be derived from the numbers given in [59]. The intra and inter events 
fluctuations of p are modelled as growing like yfp: we thus just give a and 5p the factor 
\/\ + Npu/A, though further studies might be of value to parametrize these terms in a more 
adequate manner. Notice that the value given for 6p ignores the fluctuations in the number of 
PU events from one bunch crossing to the next, but this is beyond the accuracy of our model 
here. At high luminosity at LHC, Npjj is expected to be ~ 20, which implies p ~ 10 — 20 
GeV [591E01EI]. 

Assuming gaussian distributions for these three kinds of fluctuations, one can deduce the 
Higgs width due to the presence of UE/PU0 6MuE = 2y/{AM^) - (AM)2: 

SMuE = 2y^Aa^ + A^p'^ + p'^T.^— . (61) 

PtH 

For a gaussian peak, defining a 2a width means that we keep roughly 68% of the events around 
the average, which is in correspondence with the value / = 0.68 chosen for the perturbative 
calculation. 

We now have all the important results in hand to consider both UE/PU and FSR simul- 
taneously. 



4.3 Study of the Higgs width in presence of both UE/PU and perturbative 
radiation 

The purpose of this part is to give an estimate of how one should choose the couple of filtering 
parameters {n,ri). For that, one has to convolute the effects of UE/PU and perturbative 
radiation and compute the resulting reconstructed Higgs peak width, and then minimize it 
with respect to the filtering parameters. This is highly non trivial to do analytically and 
we leave it for future work. The simple choice made here is to say that, for a given n, the 
optimal r], denoted ijopt, is the one for which the two widths are equal. This is obviously not 
true in general, but seems reasonable to obtain an estimate (figure [9]) and to understand how 
r]opt changes when we vary and Npjj. Notice that, using this method, we have to impose 
Vopt < Vsat where r]sat is the saturation point (eq. (H3]l ). because beyond r]sat, increasing -q 
makes 5Mue larger without decreasing SMpT, thus solving the equation SMpp = dMjjE has 

^^Here again, for brevity, we define SMue to be the Higgs width in presence of UE and/or PU without 
referencing the PU dependence. Actually it serves only to distinguish the width due to UE/PU from the 
perturbative width SMpt- 
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no sense in this region. Finally, we numerically minimize 
V = Vopti'iT'), with respect to n in order to find Uopt- 



6Mpj, + 6M^^, calculated at 
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Figure 9: The Higgs width due to UE/PU and loss of perturbative radiation, combined as 
if the 2 distributions were gaussians, i.e. 6Mtot = \/SMprp + 5Mf jjp when n = 3. In this 



case, r/opt, though slightly larger, is approximately given by the intersection of the 2 curves, 
at least as long as r] is not in the saturation region. 

First, we would like to understand how rjopt evolves with respect to the physical parame- 
ters. The equality SMpT = SMjje gives an equation in L = In i: 



1-e 



cte 



-2L 



-4L 



-4L 



PUE- 



PtH 



(62) 



where the coefficients Co-, csp and cs can be easily calculated using eqs. (j49l54ll60l6ip : 



Ca{n, Npu, Rbb) - ^.Q^fTTy/nRbbX 1 + 



Npu 



csp{n, Npu, Rbb) ^ O.SirnRlf^J 1 + 



Npu 



cs(n, Npu, Rbb) - 0.26tt,/^rI (l + ^ j . 



(63) 
(64) 
(65) 



If the solution of eq. (j62]) for a given n is found to be above r]sat{n, /), then rjopt = Vsatin, f) 
in order to take the saturation of 6MpT into account. We start by solving this equation 
numerically. In figure [TU] we show rjopt as a function of pt^j and Npu for different values of 
n. As it should, r]opt increases with ptjj at fixed Npi/. Indeed, if ptjj grows at fixed 77, Rbb 
decreases and so does the effect of UE/PU, whereas the perturbative radiation is kept fixed 
(no dependence on Rbb)- Notice also, for n = 3, that the values obtained for rjopt are roughly 
consistent with the choice in [15] where we had rj = mm(0. 3 /Rbb, 1/2). The saturation comes 
into effect at relatively low pt^^ , around 400 — 500 GeV. Above this value, the total width is 
small and hadronisation corrections start to become relevant, so that the results presented 
on these plots become not very reliable. However, for pt >~ 500 GeV and 77 > rjsat, the Higgs 
width due to perturbative radiation and UE/PU vary slowly with 77 and we thus believe that 
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the precise value chosen for rj is not so important: one can take any value above rjsat without 
changing the result too much. The decrease of T]opt with Npjj seems to be weaker than one 
might have expected a priori. However, in fig. [U we can see that the negative slope of the 
perturbative width is very large, and therefore increasing the noise from PU will not change 
too much the r]opt value. 
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Figure 10: The numerical solutions (points) of eq. (I62|) shown for different values of n: 
as a function of ptj 



a 



when Npu = 0, and (b) as a function of Npu when pt^j = 200 GeV. We 
also show the corresponding approximate analytical solutions (lines) derived in eq. (f73]l . 



It would be interesting to understand analytically the evolution of r]opt with respect to 
the physical parameters ptjj and Npij. Unfortunately, eq. (|62]) cannot be easily dealt with. 
That's why we have to make an approximation: in this equation, one of the 3 terms under the 
square root may be dominant when r] = rjopt- At first sight, one would expect that at low 7]opt, 
the c^e~^^ term, which scales like r]"^, should be the largest, whereas at large Npjj, it should 
be the c|e~^^ term that is the largest one as it scales like Np^. But figure [TT] for n = 3 
reveals that the csp term surprisingly brings the largest contribution to 5Mue for physical 
values of the parameters (the same holds for other values of n). Therefore, to simplify things 
a little, one can consider eq. ([62]) and put Co- = cs = 0. However, to be more general, and to 
consider the possible situation where one of the other terms might be dominant we rewrite 
eq. (I62p in the following approximate form: 



20OC 



1-e 



-4ir/3o 



CuEPuEe-P'-Rl, 



Mh 

PtH 



where p = 1 if the term dominates and p = 2 otherwise. Moreover: 



(66) 



CuE{n,Npu) 



1.2y/TT y/n\l 1 H , if the Co- term is dominant, 

if the term is dominant, 
if the cs term is dominant. 



1.67rnY'l + ^ 
0.527ri/H fl + 



(67) 



The subtraction procedure proposed in [59] seems to eliminate most of the fluctuations from the csp and 
ce terms, so that the remaining Ca term would be dominant in this case. 
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Figure 11: 5Mue computed aX rj = rjopt with respect to (a) pt^j when Npjj = and |(b)| A^p^f 
when ptjj = 200 GeV. On these plots is also represented the contribution to 5Mue of each 
term separately. When the UE/PU width falls below the saturation line SMfjE = ^M^at, 
then rjopt = Vsat- 



Eq. ()66p can be written in a slightly different way: 



with: 



^=ln( ), (68) 



BpT = 4tt/3oCpt , (69) 
BuE = l- 2poas In ( ^ „p ) . (70) 



Despite its simpler form, eq. (j68p for L cannot be solved analytically. Here comes the second 
approximation, which is to make a perturbative expansion: 

= + + ^ (("^-^) ) • ('^1) 

Neglecting the O ((a^-L)^) term, the resulting quadratic equation immediately implies 



Lopt = — 7ir-^ ■ (72) 



-BuE In 5^ + ^jBfjj^ ln2 ^ + SpQa^pBuEBpT 
Taking into account the saturation effect, r]opt is then given by: 



^ f e , if Lopt > - lnr]sat , /^^x 
1 Vsat , otherwise . 

We used this expression with Cue corresponding to the 6p term in eq. (j67p and p = 2 in 
order to plot the approximate solutions in figure [101 This reveals that the above relation for 
Vopt (eq. (f73|) ) works rather well, within a few %. 
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As a second step, we would like to find the optimal n, denoted Uopt- This also should 
depend on the way UE/PU and perturbative radiation are combined. However, as a simple 
approximation, one can combine them as if they were both gaussian distributions. Therefore, 
one should minimize 

SMtot (n) = ^5Ml^{n)+5M^j,{n) , (74) 
computed at r] = r]opt{n) for a given pt^j and Npfj. 
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Figure 12: SMtot computed at -q 
values of n. 

The results are plotted in figure [T2j We can notice that the larger n, the narrower the 
peak, and thus the better the result. However, one should keep in mind that when n increases, 
the optimal Rfnt = rjRhb becomes small, and we have to deal with hadronisation corrections 
that grow as l/Rmt |56j as well as detector resolution and granularity 6r]x6(j) = 0.1 xO.l that 
both start to have an important impact on the reconstructed Higgs width, and thus degrade 
the results presented here. In section 14.51 we will examine what happens when we include a 
very rough estimate for hadronisation corrections. However, at first sight, it seems that one 
should definitely not take n = 2. The value n = 3 chosen in [15] is good, but it may be 
possible to do better with n = 4. Beyond this value, the optimal Rmt falls below ~ 0.2 (cf 
figure [T0|) . which is too small for this study to be fully reliable, as we shall see in section ITSl 



4.4 Variations of the results with z and / 

Until now, we have only presented some results for / = 0.68 and z = 1/2, z being defined as 



z = mm 



Eh Eh J 



(75) 



with Ei the energy of particle i in the Higgs splitting into bh. What happens if we change 
these values? 

Let us start with z. Though the Higgs splitting into bh is more often symmetric than in 
QCD events (and this is what was used in [15] to distinguish it from pure QCD splittings). 
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it still has a distribution in z that is uniform in the range: 



1 



1 

1 - 



<z<\. (76) 



which follows from simple kinematics in the limit m;, = 0. But in order to reduce the large 
QCD background, one usually cuts on small z, so that 

Zcut < z (77) 

with Zcut ~ 0.1. As an example, assume the b quark carries the fraction z of the Higgs 
splitting. In such a case, b and b are not equidistant from the Higgs direction: they are 
respectively at an angular distance (1 — z)Rbb and zRbb from H (see for instance figure [16] 
in appendix |A|) . Therefore, as UE/PU particles tend to cluster around the perturbative 
radiation, eq. ([^7j) has to be modified: 

6gH ~ zRbb or OgH ~ (1 - z)Rbb , (78) 



for a given UE/PU particle g in the filtered jet. This leads to the modification of eq. (^6 
according to g being relatively close to b (region called "Ji") or b (region called "J2"): 




2p: 



tH 



\ 96^1 96J2 / 



^ -pA, + ) . (79) 



In this calculation we used eq. (fl5]) . To compute the dependence of the fluctuations on z, 
we take the simplest case n = 2. For the a and S fluctuationSjthe terms pAi and pA2 vary 
independently, leading to the following contribution to SMuel^ 



(80) 



where 

with Co- and cs given by eqs. ()63|65p for n = 1. Concerning the 5p fluctuations, the 2 terms 
pAi and pA2 vary the same way from one event to another. Therefore, if it were only for the 
Sp term, we would write pAi = pA2 leading to: 



the factor of 4 comes from the fact that we compute the width at 2a. 
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where 



with csp given by eq. (f64ll for n = 1. Adding all these contributions, 



',5p 



(83) 
(84) 



this apparently leads to an enhancement of the width by a factor of 1/z. But we have to take 
into account that the coefficients csp and cs also contain a factor i?^^ (eqs. ()63|64p ) leading to 
another factor 1/z, and thus an enhancement l/z"^ at small Therefore, we can conclude 
that the effect of z 7^ 1/2 is to broaden the reconstructed Higgs peak. Such a factor may 
partly explain the width of ~ 14 GeV that was observed in [62j, to be compared with the 
various widths found in the previous subsection (see for instance figure [T2]l . and should also 
lead to decreasing rjopt- This is illustrated in fig. [131 which was obtained with the results 
derived in appendix |Bl where we carry out the above analysis for a general n. 
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Figure 13: (a) rjopt as a function of ptjj when / = 0.68 and z = 0.2 for different values of n. 
The points correspond to the numerical determination of ijopt, found solving eq. (j62p whose 
z dependence is derived in appendix [HI whereas the curves correspond to its approximate 

analytical solutions. |(b)| JMto^ = 
for / = 0.68 and z = 0.2. 



6Mprp + 5M^^ computed at r] 



rjopt as a function of pt^ 



Now, we turn to the / value. As we explained in section [^T| the choice / = 0.68 was made 
to correspond to a 2a gaussian width, as we did for SMjje, which is somewhat arbitrary. We 
would like to estimate how the results change when / is modified. We thus also consider a 
range of values for / between 0.5 and 0.8. In this case the CpTin, f) constants caracterizing 
6MpT are changed (see for instance eq. (jlO]) ). and 6Mue is also changed, i.e. eqs. ()61|62l67p 
have to be slightly modified: 



2\/2erf-i(/) cle-^^ + c^^e'^^ 



+ c|e- 



PtH 



(85) 
(86) 



*This is valid when pt^ > 



M„ 



with Ro the radius of the initial clustering of the event, in order 



for the b and b to be clustered together. For lower pt^j , there is a kinematic cut on z and the enhancement is 
less strong. 
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where erf(x) is the usual error function: 

erf(x) : 



(87) 



Notice that the constants c„, csp and cs are left unchanged with this convention. However 
Cue becomes: 



CuE{nJ,Npu) = < 



2^/2 erf-i (/) 1 + 



2^/2erf~\/)0.87rnJl + 



Npu 



if the Co- term is dominant, 
if the csp term is dominant, 



2^/2 erf"^(/) 0.26TT^/n (l + ^] , if the cs term is dominant. 



The bands corresponding to the uncertainties on rjopt that we obtain including these modifi- 
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Figure 14: Uncertainty on rjopt when / varies from 0.5 to 0.8, for different values of n as 
a function of (a) ptjj when Npu = and |(b)| Npu when ptj^ = 200 GeV. The results for 
/ = 0.68 are also plotted as a reference. 

cations are presented in figure [T^ The uncertainty that we get, ~ 20 — 30%, is not larger than 
the precision of the whole study of this paper, which limits itself to a large- A'^c leading-log 
calculation. Notice that the variation with Npjj remains small. 

One finally observes that r]sat{n, f) is almost independent of / for n = 3. In appendix iBl 
we will show that it can be approximately written as: 



Vsat ^ e [ 1 + 0.044 ( / - i") + O n / - i 



(89) 



Because of the small coefficient of its first order correction, r]sat = e~'^'^^ is a good approxi- 
mation within less than 1% for a large range of / values. But this seems to be a coincidence 
with no deep physical reason. 



4.5 Hadronisation corrections 

It is difficult to calculate what happens during the process of hadronisation, though some 
analytical results can be found concerning jet studies for instance [56l ESl El] • In particular, it 
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was shown in [56j that such non-perturbative corrections lead to a pt shift for QCD jets equals 
on average ~ — A/i?fiitCj where A = 0.4 GeV and Ci = Cp oi Ca depending on whether it is 
a quark jet or a gluon jet. This can be translated in our study by the following averaged pt 
shift for the filtered jet: 



t/ had 



(2Cf + {n- 2)Ca) 



A 



in 



(90) 



where the second equality holds in the large Nc limit. Unfortunately, there is no result 
concerning the dispersion of the pt distribution, which is the relevant quantity to compute in 
our case. Therefore, we are going to assume that the spread is of the same order of magnitude 
as the shift. This is in principle a crude approximation, but the only aim here is to illustrate 
the consequences of including hadronisation corrections in order to emphasize the fact that 
taking n too large is certainly not a good choice. Therefore, we use eq. (j46p to estimate very 
roughly the hadronisation corrections to the reconstructed Higgs mass peak width: 



SM, 



had 



{n - l)NcA Mh _{n- l)iVcA 



Rk 



2r] 



(91) 



when z = 1/2. As before, one should know how to combine perturbative radiation with 
UE/PU and hadronisation corrections in order to minimize the resulting combined width. 
However, we simply choose to minimize the quantity 



with respect to r] and plot the resulting minimal 6Mtot for different values of n (fig. [T5]l . 
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Figure 15: SMtot including hadronisation corrections computed at rj 
(a) ptjj and |(b)| A^pj7 for different values of n. 



r]opt as a function of 



The first thing one can notice on these plots is that increasing n also increases the hadro- 
nisation corrections. For n = 5 they become so important that it is now clearly not an 
optimal filtering parameter contrary to what could be deduced from figure [T2l The relevant 
Pt region in our study is roughly 200 — 400 GeV where we find the major part of the Higgs 



27 



cross-section above 200 GeV and where our results are more reliable (see section 14. 3p . In this 
region n = 3 gives the best result. And at high PU, n = 3 and n = 4 both seem optimal, 
whereas n = 2 is far from being as good. 

To conclude, our estimates seem to indicate within the accuracy of our calculations that 
n = 2 is not a good choice, nor is n > 5. Taking n = 3 or n = 4 gives equally good quality 
to the mass peak. Increasing the hadronisation effects with respect to eq. (j9ip . would lead 
to 

n-opt — 3, whereas if we lower them, we would find Ti-opt — 4. The only thing we can say is 
that n = 3 and n = 4 both seem to work rather well. 

One way to go beyond these results would be to use event generators like Herwig [65l |66] 
or Pythia [67], to compute directly the Higgs width in presence of UE/PU, perturbative 
radiation, ISR and hadronisation, and to find for which value of the couple (n, rj) the re- 
constructed Higgs mass peak width 6M}{ becomes minimal (that would still depend on pt^ 
and the level of UE/PU). But our study here aimed to understand as much as possible the 
physical aspects behind such an optimisation, the price to pay being larger uncertainties on 
the result because of the necessary simplifications that were made. 



5 Conclusions 

This work has investigated the effect of QCD radiation on the reconstruction of hadronically 
decaying boosted heavy particles, motivated in part by the proposal of [15] to use a boosted 
search channel for the H ^ bb decay. Though this article took the Higgs boson as an example, 
all the results presented here can be applied to the W and Z bosons, as well as any new 
colorless resonance decaying hadronically that might be observed at the LHC. The main effect 
of the QCD radiation is to distort and spread out the boosted heavy resonance shape well 
beyond the intrinsic width of the resonance. The aim therefore is to calculate the resulting 
resonance lineshape. This is a function of the parameters of the reconstruction method, 
notably of the "filtering" procedure, which aims to limit contamination from underlying event 
and pile- up, but which causes more perturbative radiation to be lost than would otherwise 
be the case. 

Calculations were performed in a leading (single) logarithmic and leading colour approxi- 
mation, which is the state of the art for this kind of problem. Analytic results were provided 
up to In^ for n = 2, and all-orders analytic results for the cases n = 2 and n = 3 were 
given for the terms that dominate in the small r] limit. Numerical fixed-order results up to 

In^ and all-orders resummed results were also given and are treated in more details in 
appendix [C] for a range of n and r]. For the n = 2 and n = 3 cases there is quite acceptable 
agreement between the small-ry analytic results and the full numerical results, even for values 
of ~ 0.5. 

One unexpected feature that was observed was the behaviour of order-by-order expansion 
as compared to the resummed result: indeed there are indications that the series in Og In -^jq 
has a radius of convergence that is zero (but in a way that is unrelated to the renormalon 
divergence of the perturbative QCD series). This seems to be a general feature of the non- 
global logarithm series. Its practical impact seems to be greater for large r], or equivalently 
when the coefficients of the "primary" logarithms are small. 

With these results in hand, it was then possible to examine how the perturbative width 
of the resonance peak depends on the parameters of the filtering. Though this was accessible 



only numerically for the full range of filtering parameters, figure 8(a) lends itself to a simple 



28 



parametrisation for practically interesting parameter-ranges. 

This parametrisation was then used in section 14.31 together with a parametrisation for 
the effect of UE and PU, so as to examine how to minimize the overall resonance width as 
a function of the filtering parameters and of the physical parameters of the problem such 
as the resonance pt and the level of UE/PU. The approximations used might be described 
as overly simple, yet they do suggest interesting relations between optimal choices of the 
filtering parameters and the physical parameters of the problem. Though it is beyond the 
scope of this article to test these relations in full Monte Carlo simulation, we believe that 
investigation of their applicability in realistic conditions would be an interesting subject for 
future work. It should also be noticed that the methods used in this paper may be adapted to 
other reconstruction procedures like jet pruning [T3] and jet trimming [19] as well as filtering 
as applied to jets without explicit substructure [T8]. 
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Appendices 

A Analytical considerations on the non-global structure of the 
perturbative expansion 

We derive here all the results presented in section [2l The Higgs boson is taken to move 
along the x direction and the angular coordinates {9, (j)) are defined with respect to the Higgs 
direction (so that 9 = corresponds to the Higgs direction and </) = to the y axis for 
instance). As it is very boosted, the b and b resulting from its decay should be close to 
the Higgs, i.e. with ^(,5 ^ 1. Due to angular ordering, so will be the major part of the 
perturbative radiation from bb. To take this property into account, the angular coordinates 
{9, (p) are slightly changed into a two-dimensional vector a: 

a = 9 {cos (p, sin (/>) , 
(fa = 9d9d(p . (93) 

This vector is useful because of the property that, at small 9, one can express the angle 9ij 
between 2 particles i and j as: 

9ij = \ai- aj\. (94) 

RjjI is defined to be the vector from 6 to 6 in the a plane (see fig. [T6|) . Even if \a\ < 0(1), 
the integrations over ax and Oy are extended to 00. By doing so, the error made on a 0(1) 
result is of order O (-R^b) • 
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/"^ nRbbX 

















Figure 16: The a. plane, with various variables used in the calculation and all along this 
study. In this figure, the h quark is supposed to carry a fraction z of the Higgs energy, and 
the center of the frame coincides with the direction of the Higgs boson momentum. 



A.l Primary coefficients 

Let us go back to the integral of eq. Q and write it in the a plane: 



I{AM) 



|M(A;)pe(AM(A;) - AM)) 

Jo h 



Uhi (27r)32|fc| 

27r2 ./n k 



Q - ziJfcfcpla + (1 - z)i?fefeP 



G {AM{k) - AM) . 



For AM{k) = Mh - M one easily finds 



AM(k) = kt—A{d) , 

PtH 



with 



A (a) 



, 1 



2M 



1 

+ 1^ 



2z{l - z)Rl 2 

where we used eq. ()45|) . With this expression, the integration over kt is straightforward: 

fP'H dkt 

10 h 



(95) 
(96) 

(97) 
(98) 



P'H dkt 

MfjA(S) 



ln^+lnA(a) 



(99) 
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Notice that \a\ ~ (1 — z)Rbb or \a\ ~ zRbb, depending on whether the perturbative gluon 
emission is relatively close to 6 or 6 (due to the colHnear divergence of QCD). Thus, given that 
a Higgs spHtting is most of the time roughly symmetric {z ~ 1/2), this leads to A (a) = 0(1). 
Therefore: 

ln^(a)«ln^, (100) 
and the In^(a) term can be neglected in a leading-log calculation. One thus obtains: 



J(7?) , (101) 



with J{rj) the remaining angular integral. Introducing the coordinates (r, ip) defined in fig. [161 
J{r]) can be rewritten in the following formj^ 

m - M,) ^ i * r"'t ' W^^^T^r 4 < - < 1 . 

arccos(— T -\- T 55 COS -p 

=j„fivi r,^_ ^Jk ^ ,f„<i. (102, 

Performing the ip integration and the r one (for 7] < 1/2), one arrives at the formulae ()llll2p . 
where all the Rbb dependence is cancelled. 

One remark: if r] > 1, the b and b quarks cluster together, and the result then depends 
on the z fraction of the splitting. For instance, if r/ > 2, J{r]) can now be written!^ 

A. 2 Non-Global coefficients 

The starting point here is eq. (j20p for I^*^: 

I^^{AM) = - j dkrdk2{A-nasfQ{Ei - E2)@{ki £ Jbb)^{k2 i Jbb)W2Q{AM{k2) - AM) . 

In the same way as the primary case, using the a plane and integrating over the energies of 
the 2 gluons, one arrives at: 

1 /„, /!l,f__\\2 

tNGi 



I-^{AM) = \cpCa log (^)) S2{n) , (104) 



^^A simple shift of the a coordinates gets rid of the z dependence, for instance 5' = a — (2 — i) _Rj{, 
^"the intermediate case 1 < rj < 2 has a more complicated phase space integration and is not presented 



here. 
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with 



^bb 



+ 



+ 



\ai - zRbb\'^\\di + (1 - z)Rbb\'^\d2 - zRbb\'^\\a2 + (1 - z)Rbb\'' 
1 



\di + (1 - z)Rbb\^\ai - a2p|a2 - zRbbl'^ 



\di - zRbb\'^\ai - a2p|a2 + (1 - z)Rbb\'^ 



(105) 



In all this part, ry < 1/2 is assumed. To deal with this integral, the frame is centered around b 
for instance, and 2 quantities are computed: Stot where gluon 1 is in the jet region around 
b and gluon 2 covers the whole phase space except J5, from which is subtracted Sint where 
gluon 2 covers Jb, the jet region around b (fig. [T7|) . Therefore: 






(a) Stot 



(b) S„ 



Figure 17: regions of integration (dashed region) for gluon 2 when computing Stot (left) and 
Sint (right). 



82(1]) = 2{Stot{v) - Sintiv)) , 
where the factor 2 is for the symmetric case (gluon 1 in Jb). 



(106) 
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A. 2.1 Calculation of Stot 

Using the variables = :;^^)^^ (fig- fT6|) . Stot can be written 

Stot{'n)=--^ / Uidui / U2dU2 / #1 / #2(^ 

1 



+ 



ufn|(l + 2r/ui cos V'l + )(1 + 2r]U2 cos '(/'2 + ^^^i) 

1 

72~( 



(1 + 2ryii2 cos ^^2 + v'^uDiul - 2uiU2 cos(V'i - ^2) + ""i) 

1 

^2(1 + 2r/ui cos ipi + rj'^u1){u1 — 2uiU2 cos{ipi — 1P2) + ""2) 
Doing the angular integrations, one arrives at: 

'"^ Ui fv du2 1 

/o "1 -?7^Wi Ji U2 ul -uf 



Stotiv) = -4 / du.r^^, / -^-2— -2 ■ (108) 



2 remarks about this result: 

1. The first and second terms have divergences in ni = and M2 = so respectively when 
gluon 1 is collinear to b and gluon 2 is collinear to 6, but they cancel when adding these 
terms. 

2. The part of the integral corresponding to U2 > ^, i.e. r2 > Rhh^ is null, which can be 
simply interpreted as a manifestation of the angular ordering. 

Performing this integration with Maple for instance gives the following result: 
Stot{r]) = — - - ln2(l - r/) + ln(27?) ln(l - r/) - In 2 ln(l + r/) - 2K f dilog ^ ^ 



2 2' " ^ " ^ " ^ " y 

2Li2(l - r?) - Li2 (^) + Li2 (-r?) + Li2 (i^) - Li2 (^) , (109) 



where 



dilog(x) = ^ dt^^. (110) 

Li2(x) = \\t^^^^. (Ill) 

This is rather complicated, but that expression can be greatly simplified using the relations: 

Vx > , dilog(2;) =Li2(l -x), (112) 

dilog(-x) = ^ - 1 ln2(l + x) - Li2 {-^—\ - m ln(l + x) , (113) 
3 2 \1 + X/ 

Vx with < X < 1 , Li2(x) + Li2(l - x) = ln(x) ln(l - x) , (114) 

6 

Li2(l - x) + Li2(l - 1/x) = In^ X . (115) 
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The final answer is then: ^ 

Statin) = -\- (116) 

The remarkable point to notice is of course that Stot does not depend on ry. But no simple 
explanation was found to interpret this result. 

A. 2. 2 Calculation of Smt 

Sint must now be subtracted from Stot- The computation is similar to the previous one and 
is not detailed here. However, contrary to Stot^ a simple analytical result was not obtained, 
only the following expansion: 

4 fi 31 o ,n 1921 19 889 20589 7643 on, 

= _2r,4 _ 6r?6 - -rj' - 407?!° - ^r/^^ _ _ ——r,^^ - -—r^^^ + O(r?20) ^ 

(118) 

Using eq. ()106p with the 2 previous results, one arrives at the expression eq. ^ 



B Analytical considerations on the dependence of the results 
on z and / 

B.l Dependence on z 

It is interesting to understand analytically how rj^p^ evolves when the decay of the Higgs boson 
occurs with a z fraction different from 1/2. In fact, we obtain the same result as eqs. (|72|73p 
up to a modification of the constant Cjje (more generally, we have to modify the coefficients 
Co-, csp and cs, see eqs. ()62ll65p ). The starting point is eq. (j79]) generalized to any value of n: 

Mm -J^ 

AM ^ -^y2ai{z)pAi, (119) 

where ai{z) is either ^ or depending on whether subjet i is in the Ji region (around the 
b quark) or in the J2 region (around the b quark). This is because the UE/PU particles tend 
to cluster around the perturbative radiation, which itself is emitted close to b and b. We call 
a "configuration" the set of all the coefficients ai(z). 

The result on the fluctuations depends on which ones are considered. Let us start with 
the fluctuations originating from the a and S terms. In this case the pAi terms vary inde- 
pendently, thus leading to a contribution to SMjje similar to that of eq. ([80]) for n = 2: 

^M^...,..K| = 4 (^) ' E a.iz)Hl^, , (120) 

with (see eq. (fST]) ): 

r2 2 -2L I 2 -4L 

= • (121) 
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The coefficients and cs are still computed for n = 1 in this formula. But eq. ()120p is 
only valid for a given configuration {aj}. We thus have to average over all the 2"~^ possible 
configurations (the b and b subjets are fixed). As the perturbative radiation pattern does 
not depend on z (for z not too small), each configuration arises with the same probability. 
Therefore, if k is the number of subjets in the Ji region and n — 2 — k the number of subjets 
in the J2 region apart from the b and b subjets, we obtain: 



\ PtH / 8 z\l - zY 

We can follow a similar reasoning for the 5p fluctuations, except that the pAi terms in 
eq. (|119p vary the same way from one event to the next. Therefore, for a given configuration 
{oj}, one can deduce the following contribution to SMue'- 



=4 (^|5^)' (ta,i.}\\l, (124) 



with 5sp given by eq. (j83p . As before, we have to average this result over all the 2" ^ possible 
configurations, leading to: 



f MhPue\'^ 1 Z^"- ~ 2\ f k + 1 n — 1 — k^ 



V 2pt„ J 2"-^ ^V/c/Vz 1-2 



1 v-:^/n-2\/fe + l n-l-A;\%o 



/ Mjjp[/j; y + (n - 2)(1 - 2^)^ 

V y 16^2(1-^)2 "'p- 



Si. (126) 



One can absorb all the dependence of the resulting 6M^^ in n and z into the coefficients Co-, 
and cs, and define new coefficients c' such that 



6MuE = 2. /c^2 2 + ^2 4 + ,/2 4^^^^lPt^ ^ (127) 
V op 2. ^ 



with 



cUn,Npu,Rbb,z) ~ 0.6V^V^^^^^^^i?,,y^7^^, (128) 

, , , A/n2 + (n-2)(l-2z)2 „ / iVpr; , ^ 

c^p(n, iVp^, Rbb, z) ^ O.Stt^ t(T^) ^^'^ V ^ + ^ ' ^^^^^ 

4(n, iVpt;, i?,,, z) :^ 0.267rV^^^^Slz3i?2^ (^1 + . (130) 

With these results in hand, we can easily generalize eq. (172p for any value z of the Higgs 



splitting. One only has to modify the value of Cue- For instance, the curves in fig. 13(a) 
were obtained using eqs. (pOj - lTS]) with: 



r ( AT ^ , , v/n2 + (n-2)(l-2z)2 / Npu 

CuE{n,Npu,z) = 1.6tt^ ^^^^ _ V ^ + ' ^^^^^ 
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where we just use the dominant c^^ term. 

There is also another source of fluctuations that we haven't accounted for so far. Indeed, 
even if the values of p and of the jets area were constant, we should consider the fact that 
the filtered subjets can be either in the Ji or in the J2 region. The calculation of its effect is 
similar to that of the fluctuations in a, 6p and S. Its contribution on 5M'^^ can be written: 

5Mai.,othe. = 4(^^^)^?,he.^^ (132) 

with 

/ / N y/n^\l-2z\ 2 f Npu\ 

Cother(«> Npu, Rbb, z) = vr — ^^^^ _ — R^^ (1 + I . (133) 

But we checked that its effect on SM"^^ is negligible compared to that from the dominant 5p 
fluctuation^ and therefore we did not include it. 

B.2 Comments on the uncertainty due to the choice of / 

Let us return to the observation of section H3] that 'r]sati'n, f) is almost independent of / for 
n = 3. To understand why, we have to estimate tsatif) and Cpt(3, /) (cf eq. (|13]) ). tsat can 
be deduced from the equation: 

S(")(r? = l,t,,t) = /, (134) 
for any value of n, as tsat only depends on / (see for instance figs. ([7 ]8(a) ). To be simple, we 



can use the function T,^'^^ which was widely studied in this paper. Unfortunately we cannot 
take the primary emission result eq. (jlOp because the non-global part becomes important 
when r] = 1. However, as an approximation, we can numerically compute the second order 
coefficient 02 — —3 of S(^)(r/ = l,t) and solve: 

1 - J{l)N,tsat + aatL - / , (135) 

which simply leads to 



, J(l)A^c- VJ(l)^iV|-4a2(l-/) 

tsatif) - • (136) 

This expression can be shown numerically to give tgatif) with a precision better than 1% for 
/ G [0.5,0.8]. CpT{<i,f) is harder to evaluate. One must solve 

^^^\L,t) = f, (137) 

in the limit of large L. Using eq. (j27p . which is valid in this limit, and defining the function 
h such that h{4NcLt) = S(3)(L,t) - /, we can use Newton's method with only one iteration 
on h, evaluated numerically for ANcLt = ij^l to solve approximately eq. ()137p : 

Ml) 
h'ii) ' 

2.543 - 2.334/ , 

m,CpT{3,f), (138) 



ANrLt ~ 1 



^^This efTect is strictly null when z = 1/2. When z — 0.2, we obtain: 

SMuE,Sp ^ 16.6 (JAf^^j other ' 

for pt„ = 200 GeV and Npu = 0. 

^^Starting from iNcLt ~ 1 quickly converges to the true solution. 
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which gives Cpxi'^^f) within a few %. Therefore, eq. (fl3l) leads to 



f?sat(3,/) = e 



-0.58 



l + 0.044(/--j+O (/--) 1 1 , (139) 



which is eq. 

To be complete, the same analysis can be done when n = 2, using eq. (P0|) for Cpt{2, /), 
to obtain r/sat(2, /). We will just mention that 

limr/,„t(2,/) = e-^ - 0.85. (140) 



C Convergence of the non-global series 



In this appendix, we go back to the study of the convergence of the non-global series that 
was already briefly examined in section 13.21 trying to understand a little more what may be 
behind the observed behaviour. 



C.l Case riflit = 2 

First, we start by studying the convergence of the non-global perturbative series when n = 2 
for 7] small, where the primary coefficients are known to be enhanced with respect to the 
purely non-global ones because of the presence of large collinear logarithms (cf section [2^2]) . 
Figure [18] compares the fixed-order results to the all-orders one up to . On these 2 plots 
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t t 

(a) (b) 



Figure 18: Comparison between fixed-order (FO) expansion and all-orders result when n 
for[(^?7 = 0.1 and [(b)] r/ = 0.3. 



one can notice that the series seems to converge, as was already shown in section [321 In other 
respects, the convergence looks better when r/ is larger. This may be understood using the 
following simple explanation: if we make an expansion up to order k, then the series starts 
to diverge from the exact result when the term of order k + 1 becomes roughly of the same 
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size as the function itself. In the 'S^'^\L,t) case, using the analytic expression (eq. (I25p ). this 
can be translated to 



{k + l)\ 

with L = In i. At large the solution gives 

k + 1 
t ~ . 

So the convergence is better when k and ij increase. 



(141) 
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Figure 19: Comparison between fixed-order (FO) expansion and all-orders result when n = 2 



for (a) r] = 0.6 and|(b)|?7 = 0.9. 



Let us go to larger r], where there is no large collinear logarithm anymore, and check what 
happens. This is done for rj = 0.6 and tj = 0.9 on figure [T9l One striking feature of these 
plots is that the convergence seems acceptable up to the third order, but the 4*'^ order does 
not give as good a result, and the situation becomes worse as rj increases. Said another way, 
the perturbative expansion can be trusted until the third order, but then it starts to diverge. 
The fact that the convergence looks better for small rj may come from the dominant behavior 
of the primary series, which converges very well to a nice exponential function. However, if 
one could go to sufficiently high orders, it might be possible to observe the same divergence as 
in figure [T9l when the purely non-global coefficients become of the same order of magnitude 
as the primary ones. 

To get an idea of these coefficients, the series of the plots are explicitly written below: 



S(2)(r, = 0.1, t) 
S(2)(r; = 0.3, t) 
S(2)(r, = 0.6,t) 
S(2)(r/ = 0.9, t) 



1 - 27.57t + 368.8^2 - 3195*^ + 20200t^ - 99300t^ + O {f) , (143) 



13.88t + 87.49*2 - 334*^ + 860*"^ - 1500t^ + O {t^) 
4.656t + 6.53*2 _ g ^^s _ ^q^4 _^ ^ ^^5^ ^ 

2.320* - 2.27*2 + 1.9*^ - 12*^ + O (*^) . 



(144) 
(145) 
(146) 



For r] small, the growth of these coefficients comes from the powers of the large collinear 
logarithm L. For r] near 1, the small growth that one can start to observe at the 4*^^ order 



the derivation was done at large k, but the result seems to be reasonable even for fc = 2. 
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essentially comes from the purely non-global part. Indeed, at the 4*^ order, the primary 
coefficient is positive, and it is not the case for r/ > 0.6. This non-global growth will be 
confirmed at higher orders for the slice observable in section IC.3[ 

Let us now see what happens if the perturbative series is exponentiated. This means that 

k 

instead of plotting g{t) = 1 + ait^, which is the perturbative series, one plots e-^^*^ where 

i=l 



E 

i=l 



af and e^W = ^(t) + q (t^+i 



so that 



ci = ai 



C2 = a2 - — , 



(147) 

(148) 
(149) 



Notice that the exponentiated first order corresponds to the analytical estimate for n = 2. 
One can observe on figure[20]a nice convergence for i] = 0.1 until the 4*^* order. Concerning the 
order 5, it seems that it diverges a little, but if the coefficients of the series are varied within 
their respective errors, then it can coincide with the all-orders curve. However, one can guess 
that it should not converge at the end because the all-orders function is not strictly a simple 
exponential. This will be confirmed later with the slice. When r] = 0.9, the exponentiated 
^th Qpjjgj. surprisingly almost fits the all-orders result. Even more surprising, if the 5*'^ order 
coefficient (not shown here) is varied within its error band, it seems that the exponentiated 
fit is still improved. Is it accidental? Here again, the slice example and good sense lead us 
to answer yes but we cannot be completely sure. 



0.1 



0.01 



0.001 



0.0001 





all-orders + 

expFOk=I 

exp FO k=2 

exp FO k=3 

^ expF0k=4 

"^s. exp FO k=5 - - - - 


n=2 
Tl=0.1 





0.05 



0.1 



0.15 
t 



0.2 



0.25 



0.3 



(a) 



0.5 



0.2 



0.1 





n=2 




r|=0.9 




+: 




all-orders + 


expF0k=l 

exp FO k=2 

expF0k=3 

expFOk=4 


+*. 
+•., 

-e 



0.05 



0.1 



0.15 
t 



0.2 



0.25 0.3 



(b) 



Figure 20: Comparison between exponentiated fixed-order (exp FO) expansion and all-orders 
result when n = 2 for (a) rj = 0.1 and|(b)|ry = 0.9. 



C.2 Case Uf^t = 3 

For small enough ry, it was found analytically in section [231 that the all-orders result for the 
Higgs AM distribution can be expressed as a simple function of t (eq. (f27ll ). which is of the 
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Figure 21: Comparison between fixed-order expansion and all-orders result when n = 3 for 



(a) ?7 = 0.3 andp^ry = 



form g{4NcLt) with: 



9it) 



1 + 



* 1 
dt'- 



t' 



(150) 



|afc|<-(2^ + 0(l)) 



(151) 



The radius of convergence of the Taylor series for g is infinite as the coefficients resulting 
from its expansion can be shown to be bounded for large k by: 

1 

As the expansion of this function converges, one would expect the same to occur for the 
curves obtained numerically. However, the expansion does not converge as fast as for the 
exponential, so that the t window for which there is a convergence should be smaller than 
what was obtained for n = 2. This is indeed what is observed in figure 21(a) for rj = 0.3, if 



it is compared for instance with the plot 18(b) 



In other respects, the curves on the plot 21(b) for rj near 1 behave similarly as for n = 2, 
i.e. a perturbative series that converges until the 3^^^ order, and that starts to diverge from 
the 4*'^ order. Notice that the same behaviour is observed for n > 3. 



C.3 The Slice case 

The Slice observable, studied for instance in [H], gives an interesting example of the strange 
behaviour of the non-global leading-log series. It is simply defined by the sum of all the 
particles' energy flowing into a region Q of the phase space corresponding to a rapidity 
y G [—yo,yo], with yo a parameter of the observable (figure [22|) . Here, we work in the qq 
center of mass frame and the quarks are assumed to move along the z axis. This observable 
is non-global as shown in fig. [22l and it is interesting in 2 ways: 

• A simple change of fram^^ shows that it is more or less equivalent in the leading-log 
approximation to the AM observable in the filtering analysis using anti-Zct with n = 2 
and r] = e"^" when 77 <C 1 , while being faster to compute numerically due to the absence 
of clustering. 



'or more simply comparison between primary results, eq. (2.5) from [JT] and eq. (|lUp from this paper. 
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Figure 22: The region between the 2 dashed hnes, and the gluon's configuration leading 
to the appearance of non-global logarithms. 



• When e~^" ~ 1, which should very approximately correspond to 77 = 0(1), the strange 
behaviour of the non-global series, observed in section [3l is clearly confirmed with the 
addition of the 5^^ and 6^^ orders. 

The series are represented in figure [23] for 2 different values of i/q. The plot for yo = 2.3 is 



here for a comparison with figure 18(a), i.e. when n = 2 and rj = 0.1. One can notice that 
until order 5 the behaviours of the 2 series are very similar. The slight difference comes from 
the fact that the C/A algorithm was used there instead of anti-kt- This can also be seen 
on the expansions eqs. ()143ll52p . However, a remarkable effect is the G*'' order curve which 
does not improve the fit with the all-orders result anymore (it even makes it worse). When 



yo = 0.5 (fig. 23(b)[ ), this effect is enhanced: indeed, one notices that the 3'''^ order gives the 



best result, with the 2"^^ order even worse than the 1^* one. And here again, adding more 
orders shifts point of disagreement to smaller t. 

To get an idea of the coefficients in this case, the series of the plots are given below: 

S(yo = 2.3, t) = l- 27 M + 351. It^ - 2673*^ + 13900t^ - 99500t^ + 2 lO^t^ + 0{t^) , (152) 
S(yo = 0.5, t) = 1 - 6i - 9.108^^ + 114*^ + 1740t^ - 68400t^ + 1.6 IQ^t^ + 0{f) . (153) 

The growth of the coefficients for yo = 2.3 until the 5*'' order essentially comes from the 
powers of ^NcUo when expanding the primary result, which is the expression of the collinear 
divergence near q and q, whereas, for yo = 0.5, it essentially takes its origin from the purely 
non-global part. There is no large enhancement due to collinear divergence which can explain 
it. As an example, we also show the function S{yo,t), which is defined as in [351141] to contain 
the purely non- global part of the result: 

where, for the slice, the primary contribution S^^^ can be written as: 

s(^)(yo,t) = e-^^^^^°*, (155) 

in the large- A'^c limit. The plots for S{y(),t) are shown in fig. [25] One observes the saturation 
property noticed in [H] which leads to very similar plots for yo = 0.5 and yo = 2.3. There is 
clearly no convergence of the perturbative series. 
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Figure 23: Comparison between all-orders and fixed-order results for the Slice when 
yo = 2.3and[(b)]yo = 0.5. 
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Figure 24: Comparison between all-orders and fixed-order results for the purely non-global 



part S{t) of the Slice when (a) yo = 2.3 and|(b)| j/p = 0.5 
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Therefore, the non-global leading-log large- A'^c series seems to behave badly at high orders. 
Does this mean that it is an asymptotic series like the Standard Model is known to be [68|'i^ 
This study cannot answer such a question but, at least, one should be aware of the strange 
behaviour of the non-global series. 
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Figure 25: Comparison between all-orders and exponentiated fixed-order results for the Slice 
when (a) yo = 2.3 and [(b)] yo = 0.5. 



To finish, let us mention the exponentiated results represented in figure [25j They show 
that the convergence observed in the case n = 2 is clearly only illusory. One point to notice 
is that the case yo = 2.3 is slightly different from its ry = 0.1 counterpart as the divergence 
starts to be visible at order 4 instead of 5. This may be because using the C/A algorithm in 
the filtered Higgs mass observable reduces the impact of the non-global logarithms [42 ^ B3] 
and the primaries still dominate at order 4, whereas it is not the case anymore for the slice. 
Notice also the very nice fit given by the exponentiated 2"°' order curve when yo = 0.5. 
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